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Abstract. In this article we deal with the approximation of solutions of non- 
linear evolution equations of the form A(u(t)) + f{u(t)) = u'{t), the numerical 
analysis of solutions to this problems will be performed using some methods 
from particular algebras of operators which are sometimes represented by uni- 
tal subalgebras of the unital C*-algebras of operators that are generated by 
some basic operators say l,a, D(-) g C{H™{G)) that in some suitable sense 
are related to the operator A{-) £ C{H"^{G)) in the evolution equations, par- 
ticular cases where the operator algebras do not verify the C*-identity with 
respect to the norm chosen are also studied, when applicable basic C*-algebra 
techniques are implemented to perform some estimates of numerical solutions 
to some types of problems, in all this work expressions like H"^(G) will repre- 
sent a prescribed discretizable Hilbert space with G CC K". 



1. Introduction 

In this work we will focus our attention in nonlinear evolution equations of the 
form: 

A{u{t)) + f{u{t))^u'{t) 

with u £ H™-{Gt) and where H™- [G) , H"^ [G t) are prescribed Hilbert spaces with 
G CC M" and := G x [0,T],r > 0. In the last equation we have that A{-) e 
C{h^(G)) is an operator variable in time that represents in some suitable sense the 
boundary conditions of the the initial boundary value problem described by (11.11) . 

We will call any C*-subalgebra ^(iJ"(G)) of the unital C*-algebra £/{H"{G)) 
of operators generated by l,a G C{H'''\g)) such that e £/(iJ™(G)), a 

particular C*-algebra. 

In this work, to obtain discrete representations of the operators involved in the 
semilinear boundary value problems, we will use an abstract procedure known as 
particular representation method, that is based in the ideas presented by Stein- 
berg and Robidoux in jl4| and |17) , and is applicable to compatible discretization 
techniques. 
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2. Background 

In this section we will present some important concepts that are needed to deduce 
many of the results presented in this work, the first concepts presented will be those 
related with sesquilinear forms of operators in prescribed Hilbert spaces iJ"(G). 

If in a prescribed Hilbert space H"-{G) we have that the inner product (•, ■)h"(g) 
is related to a map A^[-](-) : H"-{G) x H"{G) -> C by the expression: 

{^:y)H^iG)-=-^My) (2-1) 

then for any operator A e j2/(_ff"(G')) with s/{H'"-{G)) a particular C*-algebra of 
operators we can define the related sesquilinear form by. 

Definition 2.1. Sesquilinear form of an operator. For any operator A G with 
£/ :— £^{H^{G)) a particular C*-algehra of operators in a prescribed Hilbert space 
H'^{G) the form A[x]{y) defined by 

A[x]{y) := M[Ax]{y) = (Aa;,y)^„(e) (2-2) 
will be called sesquilinear form of the operator A. 

2.1. Contractions. In this section we will present some important definitions and 
results related with contractions wich are operators described by the following 

Definition 2.2. An operator T G B{X) with X a Banach space, for wich 

\\T{x)-T{y)\\ < \\x^y\\,x,yeX (2.3) 

is called a contraction. If there is a K < 1 for wich \\T{x) — T(y)|| <K\x — y\, T 
is called a strict contraction. 

Theorem 2.1. Contraction Mapping Principle. A strict contraction T £ B{X) on 
a Banach space X has a unique fixed point, ie., there exists a unique x £ X such 
that T[x) — X. 

Proof. First let us prove uniqueness. If T{x) — x,T{y) — y, then ||a:: — y\\ = 
\\T{x) - T{y)\\ < -ft: ||x - 2/11 . Since < 1 and \\x - y\\ > 0, we conclude \\x - y\\ = 
0, i.e. X — y. To prove existence, we first note that T is automatically continuous 
since ||a; — y\\ < K~^£ implies ||T(a;) — T{y)\\ < e. Now, let xq be arbitrary and let 
Xn ■= T"{xo). We will show that {x„} is Cauchy. 

||a:„ - a;„+i|| = ||r(a;„_i) - T(a;„)|| < - a;„l| (2.4) 

< \\Xn-2 - Xr^-iW (2.5) 

••• <i^"-i||xi-xo|| (2.6) 

Thus if n > m, 

n 

||a;„-a;„l| < ^ \\xj - Xj^i\\ (2.7) 

j—7n-\-l 

< K"'{1- K)-^\\xo-xi\\ ^0 (2.8) 

n— f oo 

Thus Xn is Cauchy, so a;„ — > a; for some x. Since T is continuous, T{x) = limTa;„ = 
lima;„+i = x wich proves the theorem. □ 

From expressions (|2.7p and (|2.8I) in the proof presented here, taking n oo, we 
can obtain the following. 
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Remark 2.1. A fixed point x E X of a contraction T g B{X) in a Banach space 
X , satisfies the following estimate 

\\x -x^\\< - K)-^ llxi - xoW . (2.9) 

2.2. Elliptic Operators. A differential operator defined for any prescribed u G 
H^iG) with G CC M" in the form: 

Au :— ''^di{ai,j{x)dju) + ^'^aj{x)djU, x e G (2-10) 
with di :— dxi, can be expressed in the following sesquilinear form 

A[u]{v) J {^2{aijdju)diV + '^^ajdjU -v}, x&G (2-11) 

that will be related to uniformly elliptic equations by the following. 

Definition 2.3. Strongly elliptic sesquilinear form. A sesquilinear form like (|2.1ip 
is said to be strongly elliptic if there is a constant cq > such that 

i?e^a,,,(x)^,J, > co^ ie,f ,xeG,^e C". (2.12) 

id j 

If S!/{H'^{G)){G) is a particular C*-algebra on a Hilbert space H^{G), the 
sesquilinear form £ of a uniformly elliptic operator L £ j2/(_ff^(G))(G)) can be 
expressed in the following way 

C[u]{v) ~ j OijdiU ■ djV + o,jdjU ■ v} (2-13) 
sesquilinear forms of this type are said to be 

Definition 2.4. Strongly Elliptic Form. A sesquilinear form like (|2.13p is said to 
be strongly elliptic if there is a c^ > such that 

ReJ2^^,MJ><^'LY.\^J\', x€G,i&C^. (2.14) 

Now we will introduce some important definitions related to the well posedness 
of boundary value problems and that generalize the above results to Hilbert spaces 
i/"(G) and particular C*-algebras ^(iJ"(G)) of operators. 

Definition 2.5. Coercivity of a form. The sesquilinear form A[-]{-) of an operator 
A G ^(i7"(G)) with £/{H"{G)) a particular C*-algehra over a Hilbert space H"-{G) 
is said to be H^[G)- coercive if there is a ca > such that 

\A[x]ix)\ > caMI^^c^ ,x e H-{G) (2-15) 
with \\x\\ f^Q-^ {x, x) j^^^^y 

Definition 2.6. Ellipticity of a form. The sesquilinear form A[-]{-) of an operator 
A G ^(i7"(G)) with £if{H'^{G)) a particular C*-algehra on a Hilbert space if"(G) 
is said to be H^[G)- elliptic if there is a ca > Q such that 

ReA[x\{x) > CA \\x\\\^,G) ' ^ e i?"(G). (2.16) 
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Definition 2.7. Bounded Operator. A symmetric V -elliptic operator A G .s/{V) 
with H^{G) ^ V is said to be V -bounded if there is a constant < Ka < oo such 
that |^M(^;)| < Ka \\u\\y \\v\\y, u,v eV. 

An operator A e ^{H'^{G)) is said to be TZ-elliptic witli i?"(G') ^ V ii its 
sesquilinear form is ^-elliptic, in a similar way if a sesquilinear form ^[•](-) 

is y-coercive ttien ttie related operator A e ^(i?"(G)) is said to be ^-coercive. If 
for a given Hilbert space H^{G) and a prescribed subspace H"'{G) ^ V we have 
a symmetric y-elliptic operator A S £/{H^{G)) with £/{H^{G)) a particular C*- 
algebra generated by a given operator a e J^{H'^{G)), and if we have that A — a^a, 
it is not very difficult to see that 

A[x]{y) = (2.17) 

= (a^aa;,?/)^„^^^ (2.18) 

= {ax,ay)fjr^^G) (2-19) 

- {^,y)a (2-20) 

the equation (|2.20p is often called a-inner product and is useful to perform several 
estimates related with this kind of operators. A particularly important example 
of this type of inner products is the corresponding to the laplacian operator A :— 
V • V = J2k with V —: grad — [di, - ■ • , 9„] and V- =: div in particular function 
spaces, one of this examples will be presented now. 

The importance of previous properties will be remarked in the following 

Lemma 2.1. Lax-Milgram Lemma. For a given operator A S £/{H) with H a 
prescribed Hilbert space if we have that A is V -coercive and V -bounded, with H ^ V , 
then for each given f (zV , there exists a unique u £V such that 

A[u]{v)^{f,v)y, v^V (2.21) 

Furthermore, if for a given f the solution u of (|2.2ip is denoted by 

u^gf (2.22) 

then G e L{V). 

Proof. For any u,v &V and f £ H, with V ^ H, it can be seen that ^[•](-) induces 
a scalar product in the following way 

{u,v)y -.^ A[u]{v),u,v eV (2.23) 

and also we can write f(v) :— {f^v)^ since A is ^-bounded and 

{u,v)y^ f{v),v£V (2.24) 

by the Riesz representation theorem there exists u £ V such that {u,v)y — 
f{v),yv G V, hence there exists a solution to (|2.2ip . on the other hand we have 
that if there are two such solutions ui, mi E V then 

A[ui]iv) = {f,v)H = A[u,]{v),v G V (2.25) 

wich implies that 

A[ui-U2]iv)^0,veV (2.26) 
since this is true for each v £ V taking w = ui — U2 we get 

= A[ui - W2](wi - U2) > OA \\ui - ually > (2.27) 
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wich implies that ||ui — U2\\y = hence ui — U2- Uniqueness ahows us to define a 
map Q € C{V) such that for a given / e F we have that Gf is the unique element 
of V that satisfies 

A[gf]{v) = {f,v)y,v(,V (2.28) 
clearly for /i, /2 G V and scalars ai, 02 G C we will have that for every v 

^[01^/1+02^/2]^ = aiA[gh]{v)+a2A[gf2\{v) (2.29) 

= arA[h]{v) + a2A[f2]{v) (2.30) 

= ^[al/i+a2/2](i^) (2.31) 

Thus Q{aif\ + 02/2) is solution of 

A[y]{v) = {aih + a2f2,v)y = AiaiGfi + a2gf2Kv),v G V (2.32) 

hence G G jC{V). □ 

Example 2.1. // we take A = —Ajj2(^q-^ := — A|jj2((3) with G := [—1, 1]^ and with 
a := V and := V- and if we take the scalar product {u,v)j^2(^q^ to be defined by 



{u, v)q := / uv (2.33) 
Jg 

by the first Green's formula we can obtain that 



{Au,v)i^2(G) 


(2.34) 




(2.35) 


/ Vu-Vv 


(2.36) 


Jo 






(2.37) 



now we can see that by the Poincare inequality there is a number c{G) > that 
depends on the Lebesgue measure of G such that 



= c{G){Vu,Vu)^2^G^ (2.39) 
= c{G)\A[u]{u)\ (2.40) 
= c{G)\\u\\l (2.41) 



wich implies that 



^\\m(G) = MIhg) + \\^^\\1^(g) (2-42) 
< c{G)'\\Vu\\'+\\Vu\\l.,G. (2.43) 



1(G) — ll"llL2(G) -r II V u.||^2(G) 

Il2(G) + ll^"llL2(G) 

= (l + c(G)2)||V^*||^.(c) (2-44) 

= {l + c{Gr)\\u\\l (2.45) 

= {1 + c{G)^)\A[u]{u)\ (2.46) 

wich means that there is o = (1 + (c(G))^)^^ > for all u G Hq{G) such that 
A[u]{u) > ca hence we have that A is H^{G)-coercive. On the other hand 
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using the Schwarz inequality we can check that for w G Hq (G) we will have 

\A[u]{v)\ = {au,av)^2^G) (2-47) 
< Il«w||i2(c;) ||at;||i2(G) (2.48) 



< 



i\^^\\l^{g) + M\Ihg)\ [\\"''"\\h{G) + \Mh{G)\' (2-49) 
= \MHHG)\Hm(G) (2-50) 
therefore A is Hq{G) -bounded. 

2.3. Parabolic Operators. A differential operator P 6 £^{H"{Gr)) defined by: 
Pu := Au-dtu (2.51) 
= ''^di{aij{x,t)dju) + ''^aj{x,t)djU — dtu (2.52) 

i,j 3 

with {x, t) G Gr = G X [to, T], and where A is l^-elhptic for H"-{Gr) ^ V, is said to 
be parabohc. In this section we will introduce an important condition trough the 
following. 

Definition 2.8. Dissipative Operator. An elliptic operator A e Ji'{H'^{G)) is said 
to be dissipative if we have that 

ReA[x]{x) <0 ,x e dom{A). (2.53) 

The solution to a parabolic equation can be expressed using a general concept 
presented in the following definition. 

Definition 2.9. Contractive Semigroups. A set {st :t> 0} of operators on H'^{G) 
that satisfy the following conditions 

SI: \\st{x)\\ < \\x\\, X G dom{A),t > 

S2: St+r{-) = StO Sr{-), So = l,t,T> 0, 

S3: S(.)(x) e C([0,oo)),F"(G)), x G if"(G). 
called contractive semigroup conditions. 

A semigroup {st : t > 0} is said to be generated by an elliptic operator A G 
^(iJ"(G)) if we have that 

lim h-\sh{-) - 1) = A{-) (2.54) 

i->0+ 

If we have a dissipative elliptic operator A G (iJ" (Gr)) a solution v G H"{Gt) 
to the Cauchy problem 

Av = dtv, {x, t) GGr = Gx [0, T] 

v{x,0) =vo,x gG (2.55) 
Bv = Vb,x G dG 

will satisfy 

Dt{\\v{x,t)\\) = 2Re{dMx,t),v{x,t))jjn^G.) (2.56) 

= 2ReA[v]{v) <0, t>0, (2.57) 

so it follows that ||f(a;, i)|| < ||w(a;, 0)||, t>0. This shows that 

WstVaW < \\vohvo e dom{A),t> 0, (2.58) 

so each st is a contraction in the H" (G)-norm. 

The relations presented above can be sumarized in the following result. 



1/2 



|2 



1/2 
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Theorem 2.2. If A G C{H^^{G)) where H^^[G) is a discretizable Hilhert space, 
and if A is closed, densely defined and dissipative then it generates a contractive 
semigroup. 

If we represent parabolic equations by abstract evolution equations of the form 

v'{t) ^ Av{t) + f{t) (2.59) 

with A e C{H'^ (G)) , and where the Cauchy problem consists in finding a function 
V e C{[0,T],H"{G))nC^{[0,T),H"{G)) such that, for t > 0, v{t) G dom{A) and 
(12.591) holds, and v{0) — vq, where the initial value vq is prescribed. It can be seen 
that the expression 

^tA f M! (2.60) 

k=0 

satisfies the second and third contractive semigroup conditions, if we have that 
A G C{H"'{G)) is dissipative and —A G C{H^{G)) is symmetric V^-elliptic, then it 
can be verified that we will also have 

||e*^w|| < I J2 ^-Tr~l ll^ll = e""* Ikll < t>0,ve dom(A) (2.61) 

fc>0 

with TO := infj^^o ll^^ll ^ [-^[^^K^^)!, this clearly implies that {st : St '■— e''^,t > 0} is 
a contractive semigroup. Using the last results we can obtain the following useful 
results. 

Corollary 2.1. If A is the generator of a contractive semigroup, then for each 
Uq G dom{A) there is a unique solution u G C^([0, cxd), //"(G)) of (j2.59l) with 

u{Q) = Uq. 

Theorem 2.3. If A is the generator of a contractive semigroup, the for each uq G 

dom{A) and each f G {[0, oo), II^{G)) there is a unique u G ( [0, oo), iJ"(G)) 
such that u{0) — uq, u{t) G dom{A) for t > 0, and 

u'{t) = Au{t) + f{t), t > 0. (2.62) 

Proof. It suffices to show that the function 

g{t) = / St^rf{T)dT, t > 0, 

Jo 

satisfies (|2.62|) and to note that g{0) — 0. Letting z = t — r we have 

h-\g{t + h) ~ g{t)) = f s,{f{t + h-z)~ f{t~z))h-^dz 

Jo 

i-t+h 

+h^^ / S2/(t + h - z)dz 

so it follows that g'[t) exists and 

g'{t) = I sj'{t-z)dz + stf{0). 
Jo 
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Furthermore we have 

h-\g{t + h) ~ g{t)) = h^^y^ +hst+h-r.fiT)dT- St-rf{T)dT 

ft 



{Sh - l)h-^ f St-rf{T)dT 

Jo 

rt+h 

+ h-' / St+h-rfiT)dT. (2.63) 



Since g'{t) exists and since the last term in (|2.63l) has a hmit as /i ^> 0+, it follows 
from that 

St-T f {T)dT e dom{A) 



and that g satisfies (|2.62p . □ 

Example 2.2. Given the initial boundary value problem described by 

dlu = dtu, X e (0, 1) 

ulo,t) = u{l,t) = (2.64) 
u{x,0) = uo{x),uo e H^{[0,1]) 

we can verify that 

{dlu,v)^ = / dlu-v (2.65) 

' " J[o,i] 

= dxu ■ dxv (2.66) 

J[o,i] 

= - ("'«)ifi([o,i]) (2-67) 

wich clearly implies that ReA[x]{x) < hence A = is dissipative and therefore 
{st : St — e*^== , t > 0} is a contractive semigroup over L^([0,1]). It follows from 
T \2.1\ that p.64p has a unique solution u £ Hq{G x [0, oo)). 



3. Discretization Schemes 

In this section we deal with the discretization processes applied to operators 
involved in partial differential equations wich permits us to obtain the discrete 
particular C*-algebras £/n,h{H"^{G)), with H'^iG) a discretizable Hilbert space, 
see [IH], first we will present the following definitions. 

Definition 3.1. Grid: For two given sets G C and G = {0, • • • , Af„} C Z, 
with Mm a number that depends on a fixed number m, a fixed value h G and a 
bijection f : Z : G, 3 k ^ g £ , the set G,n,h = {9k & G : gk ^ f{k),k £ 
G} is called a grid in G of size h and length Aim or simply a grid. 

For a given discretizable Hilbert space H"^{G) one can define an operator Pn^h £ 
C{H^{G)) called particular projector and defined in the following way. 
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3.1. Particular Representation of Operators. 

Definition 3.2. Particular Projector: An operator Pn^h G C{H™ (G)) , with 
H"^{G) a discretizable hilbert space, that satisfies the relations: 



Plh - Pn,h (3.1) 
Pn,hX > X (3.2) 

||^n,h-l|| <C(.)/l^" (3.3) 



will he called a perticular projector, in i3.3\) ||-|| represents any prescribed norm in 
C{H) and Hm is a number that depends on m that will be called projection order 
with respect to ||-||. 

A particular projector g C{H™{G)) can be factored in the form P„ /i = 
Pu,Hpl^f,, where G C{H^{G), (if-;,(G))*) and p,,,^ £ £((i?-JG))*, 
H™f^{G)) are cahed decomposition and expansion factors respectively. The expan- 
sion factor is related to a given grid of G by the expression p^u := {ck{u, Gn,h)} ~ u 
and the operator Pn,h is related to a basis 3^ := {pk} through Pn,hPk = Pk^Pk £ 
Using the factors described above we can represent any particular projector 
over a discretizable Hilbert space X{G) using the following diagram. 

X(G)^^X,„,„(G) (3.4) 

Pm.h 

Using particular projectors we can define the particular representation of the 
inner product map in H^"{G) and also of a given operator B £ C{H'^{G)) as 
follows. 

Definition 3.3. Inner product matrix. For a given discretizable Hilbert space H 
whose inner product is induced by the inner product map A4 € C{H,H*) in the 
following way 

{x,y)H:^M[x]{y) (3.5) 
one can define a particular representation given by 

MrnAplxKplnX) '-^ M[P,n,k-]{Prn,h-) (3.6) 

that will recive the name of inner product matrix. 

A very important property of inner product matrices will be presented in the 
following. 

Theorem 3.1. Every inner product matrix form is symmetric and positive definite 
(SPD). 
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Proof. It can be seen that for a discretizable Hilbert space H and a given particular 
projector Pm,h in H, with basis ^ = {pi, ■ ■ ■ ,Pn^}, we will have that 

{Mm,h)ij = Mm,h\pl,^hPi]ipL,hPj) 
= {PhPilH 



= Mm,h\pl,,hP3]ipL,hPi) 

and this implies that A4m,h = h- since 

we will have that Mm,h[x]{x) > for each ^ x G H\Ker Pm,h- D 

This leads to the following 

Corollary 3.1. Every inner product matrix form is invertible. 

Properties of inner product matrices permit us to express the inner product and 
corresponding norm in Hn,h for any x,y € H hy 

y^Hn^h = i'>^n,hX, mn,hy)2 = {mn,hy)* mn,hX = y* Mn,hX (3.7) 

and 

l|2^lljf„,^ = '<^^^))il,h = {l^n,hX,mn^hy)2 = \\'mm,hX\\2 (3.8) 

respectively, with m,n,h the formal square root of mathcalMn,h- 

Definition 3.4. Particular Representation of an Operator. For a given operator 
B e C{X,Y) being X,Y discretizable Hilbert spaces and being Xm,h,ym,h the sub- 
spaces relative to the particular projectors Pm,h € t^{X),Qm,h G ^(X); operator 
Bm,h e C{X^ f^,Y^ f^) given by 

Bm,h ■■= Ql,hBPm,h (3.9) 
will be called particular representation of B. 

Once we have computed the particular representation of a given operator B e 

C{X,Y) over a discretizable Hilbert spaces X,Y, in prescribed subspaccs X„i^h ^ 
X,Ym,^h ^ Y determined by a particular projectors Pm,hTQm,h, we will define the 
approximation order of a particular representation as follows. 

Definition 3.5. Approximation order of a particular representation. We say that 

the particular representation Bm,h G ^i-^m h-i h) '^f operator B £ C{X, Y) is 
of order (with ;y„, a value that depends of the prescribed number m) with respect 
to a given norm \\-\\ inY if for each x G X there exists c that does not depend on 
h such that: 

Bm,hPl,hX - ql.Bx <c,h"- (3.10) 
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Based on the definition presented above and the properties of tlie inner product 
matrices if we denote by Mm,h{X) and Mm,hiy) the particular representation of 
the inner product matrices of the discretizable Hilbert spaces X and Y respectively 
we can obtain an important definition for |H|£(x^ ^ ^) in the following way. 



\A\ 



:= sup 



sup 



\mra,h{X)v\\^ 

||mm,fc(F)Am„,/j(X)-iz| 



(3.11) 
(3.12) 
(3.13) 

= \\mmMY)Amm,hiX)~'\\2 (3-14) 
wich permits us to mimic important properties of operator algebras in C{X) wich 
is really important when we deal wich C*-algebras of operators in discretizable 
Hilbert spaces. When there is no confussion we will just refer to the norm in the 
corresponding Banach space (discrete or not) simply by As a consequence of 
the Gershgorin theorem it is not difficult to derive the following. 



sup 



Corollary 3.2. \\A\\ < 



,VAe<„(G). 



Definition 3.6. Discrete Particular C*-a,lgcbra of Operators. If a particular C*- 
algehra si is generated by finitely many operators say {l,a, 6, •••} C C(H"{G)), 
then the C*-algebra s^n,h defined by s^n,h ■= {b € C{H™^{G)) -.b := p'^ h^Pn,h,b € 
£/{H"^{G))} will be called discrete particular C*-algebra or simply particular C*- 
algebra generated by {1, an^h, brn,h, ' ' • }• 

3.2. Spectral Estimates. We will begin this section defining a very useful condi- 
tion described as follows. 

Definition 3.7. Stability. A sequence {An^h}o<h<ho! n > no of m x m matrices 
An,h is said to be h-stable if the matrices An,h are invertible for any fixed n> no € 
Z+ and all sufficiently small h, say < h < hi for any fixed n e N, and if 



sup 

m+3h<ho 



a: 



< 00. 



(3.15) 



If B is not invertible we put \\B ^ I 



the sequence {An^h}n>na.h<ha is /i-stable if and only if 

< oo, for fixed n €N. 



lim sup 
h^o+ 



00. With this convention, we can say that 

(3.16) 



Using the expressiones presented above we can obtain the following result. 

Theorem 3.2. If the sequence {An^h}n>no,o<h<ho ofmxm matrices An^h is stable, 
then A is necessarily invertible. 



Proof Let A^)^ < M for < /i < /iq and n G N fixed. Then if a; G if with H a 
discretizable Hilbert space and < h < ho, 

\\P„,hX\\ = A-^f^An,hX < M \\A„^hX\\ = M \\Pn,hAP„,hX\\ , 
\\Pn,hX\\ = \\{K,hr'K,h4 < M \\Alhx\\ = M \\Pn,hA*Pn,hX\\ , 
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and passing to the limit /i — > 0+, we get 

< Af ||At|| , ||a;|| < M\\A*x\\ 
for every x ^ H . This shows that A is invertible. 



(3.17) 

□ 



Now we wih present a useful definition and sonic important related results and 
diagrams. 

Definition 3.8. Exactly Factorizable Operator. An operator A € £/{X{G)) with 
£/{X{G)) a C*-algebra of operators over a discretizable Hilbert space X{G), that 
can be factored in the form A = a^a, with a G £/{X{G)) is said to be exactly 
factorizable in S!/{X{G)) and will satisfy the following estimates. 



Exactly factorizable operators can be represented by the following diagram. 



X{G)—^X"{G) 



X'{G) 



(3.18) 




If X{G) is a discretizable Hilbert space, a particular representation Omji G X^ hi^) 
of a G £/{X{G)) with respect to a particular projection Pm.h £ 'Pi^iG)), being 
V{X{G)) the space of all projections over X{G), can be expressed using the follow- 
ing diagram. 



x{G) ^ x:^,sG) ^ x::^,{G) ^ x:„,,(g) 



(3.19) 




Pm.h 




Xm,h{G) — X'^j^{G) 



Using diagrams p.lSp and (|3.19l) we can express the particular factorization of 
an exactly factorizable operator A G £/{X{G)) by the following diagram. 



X'*jnh{G) ^ -'^m./i(G') — -> X'*^^f^{G) 



(3.20) 




X'm.h{G) — -> X"„ f^{G) —4 X",„^(G') 

•^m,h 

Now we will present the following results applicable to exactly factorizable opera- 
tors. 



Lemma 3.1. Let A G Si/{X{G)) be exactly factorizable in s^{X{G)). Then 



\A[x]{x)\ >e\\x\Y 



Proof. It can be seen that 



\A[x]{x)\ 



I {Ax, x) I 
I (a) ax, x) 
I (ax, ax) I 
llaxf 



(3.21) 

(3.22) 
(3.23) 
(3.24) 
(3.25) 
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wich implies that |A[a;](a;)| is (SPD) and that there exists ^ > such that 

||aa;|| > (5 ||a;|| . (3.26) 
taking e = concludes the proof. □ 

Theorem 3.3. Let A G s^{X{G)) be exactly factorizable in s^{X{G)). Then we 
will have 

\\An,h\\<Uhn>l, lim = mil 



Proof. Lets define 



< |U"M| ,n > 1, lim 



\A- 



. ^ {Ax,x) ,^ {Ax^x) 
m := inf V ' / ,M := sup \ / 
0:7^0 {x,x) cc#0 {x,x) 



mn,h ■ 



:= inf 



{AP„.„x,x). 



v^O {x,x}^^f^ 
By assumption, m > £ > 0. We have 



Mn,h ■= sup 



{APn,hX, x)^^f^ 



and , analogously, 



. , {APn,hX,x)^ ^ . {Pn,hAPn.hX,Pn,hX) 

mi ; ; — = mt ■ — - 

x^O {x,X)^ f^ x=iO {PnJiX, Pn,hX) 

. r {APn,hX,Pn,hX) ^ . „ {Ax,x) 

mf , > mf — — = m 

^5^0 {Pn,hX,Pn,hX} xjtO {x,X) 



{APn,hX,x)^^ {P^^hAPn.hX,Pn,hX) 
Mn,h = sup r = sup — — ^ 

x^O \X,X)^f^ x^O \Pn,hX,Pn,hX) 

{APn,hX,Pn,hX) - {Ax,x) 



sup 



< sup ■ 



xjtO {Pn,hX,Pn,hX) xjiO \X,X} 



M. 



(3.27) 
(3.28) 

(3.29) 
(3.30) 

(3.31) 
(3.32) 



Because 



||A||=M, \\An,h\\=Mn,h, \\A-'\\=l/m, A-\ = l/m„,,., (3.33) 



we arrive at the inequalities ||yl„^;j|| < ||A|| and 
< /i < 1. It is clear that ||Ara,?i|| ||A||, and since 



< \\A i|| for all n and 



\A~^\\ < liminf 



< lim sup 

h->-0+ 



A-^ 



< A- 



it follows that 



->■ ||A i||, for each n > 1. 



(3.34) 
□ 



If we denote by Xmin{A.) and Xmax{A.) the minimal and maximal eigenvalues of 

an operator A E £/ {X) respectively, with jz/ (X) a C*-algebra of operators over 
a Banach space X, an operator A € ^(iJ"(G)) that is exactly factorizable in 
s^{H^{G)) satifies the following estimate. 

Theorem 3.4. If A e s^{H'^{G)) is exactly factorizable in s^{H'^{G)), then 

m < A„i„(A„,fc) < \^ax{K,h) < M, (3.35) 
lim \min{An,h) = m, lim A„„^(A„^j) = M (3.36) 

/i->0+ /i->0+ 
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{m,M} C a{A) C liminf ct(A„ C hmsup a{An h) C [m,M]. (3.37) 

Proof. The validity of (|335)) and (|336)) was established in the proof of T|3H The 
only nontrivial part of (j3.37l) is the inclusion 

o-(A) C liniinfCT(A„ft). (3.38) 

So assume A £ R is not in liminf cr(j4„.;i). Then there is an e > such that 

C/e(A) n cr(A„j,) = for all n > no, (3.39) 

where f7e(A) := {z e C : l^; — A| < e}. Hence f7e(0)ncr(A„^ft —Al) for all n > uq, and 
since {A ~ Al)^^ is symmetric, and therefore the norm coincides with the spectral 
radius, we get 

||(A-A1)"^|| < 1/e foraU n>no. (3.40) 

It follows that {An^h — Al}„>„o is stable, and therefore A — XI must be invertible 
by theorem T I3.2I Consequently, A ^ cr(A), wich completes the proof of (13.381) . □ 

If an exactly factorizable operator A G £/{X{G)) is X(G')-coercive and X{G) is 
discretizable we can verify that 

\AmJi[S:]{x)\ {PmJiAPm,hX,PmJiX) ^ {a^aPmjiX,Pm,hx) (3.41) 

= {aPmMX, aPrn,hx) = {Pm.haPmMX, Pm.,haPmJiX) (3.42) 

= {ara,hX,am,hx)x^^(^G)-'"^\\^\t (3.43) 

this leads us to the following. 

Theorem 3.5. // an exactly factorizable operator A e £^{X{G)) is X{G)-coercive 
and X{G) is discretizable then A^.h '■— pI^ i^Apm,h is invertible. 

Proof. From (I3.43P it is clear that Am.h is positive definite hence invertible. □ 

3.3. Discrete Time Integration. In general, evolution equations of the types 
studied here can be expressed in the form 



u'{t) = f{u{t)) 

U(0) = Mo 



(3.44) 



if we define HgiG; r) by 

Hl{G;r) -.^ {v e C{[-6,d], H\G)) : w e S^(0), Vt e [-6,6]} (3.45) 
and taking the operator T : Hg{G; r) H^{Gr) to be defined by 

T : nliG-^r) ~^ H\Gr) (3.46) 

: v^v„+ f f{v{T))dT (3.47) 

and if / G G"^^{'H'^{G; r)) then it can be seen that if we define IHl-Hfc(G r) 

\HnUG;r}-= sup ||«|| (3.48) 
te[-s.5] 

where like in the above sections ||-|| :— \\-\\h^{g)' from these expressions we get 



\\T{v)-vo\\ = 



f{v{T))dT 



< MfS (3.49) 
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where Mf :— sup^g^fc^g.^j so if MfS < R, then T{v) £ nliG;r) so that 

T e CinliG; r)). Also we have 

ft 



||T(m) -T(w)||„fc(Q. s = sup 
* te[-s,s] 



(/Kr))-/(«(r)))dT 



(3.50) 



< c/ dr (3.51) 

< c/(5||w-u||^k(G;r) (3-52) 

so if we choose 6 < mm{R/Mf,l/cf} we can see that T e C{'H^{G]r)) will be a 
strict contraction. The above results con be sumarized in the following. 

Theorem 3.6. Picard Existence theorem. Let Y{G) be a Banach space. Suppose 
that f e G'^=Hys{G;r)),ys{G;r) C Y{G). Let Mf := snp,^y^(^a;r) \\fiv)l The 
initial value problem p.44p has a unique local solution u(t). This solution is defined 
for t G {-6,5) with 5 = niin{i?/M/, 1/c/}. 

Proof. From ((3^ and from (|330)) to (|332)) we have that T G 3^5 (^; V) is a strict 
contraction for 5 < R/Mf, from this fact the result follows. □ 



Corollary 3.3. The solution u{t) to p.44p satisfies the estimate 



{cfSf 



\\u{t) ~ Ufc(t)||„.(c;r) < Y^J^fh^- (3-53) 

Proof. Follows from remark R I2.il □ 

Time integration of space discretized evolution equations of the form 

x'{t)=F{x{t)) 

x{0) = xo ^'^■^^> 

where F G C"^^(f2), 57 C can be performed using the operator S : x M — s- 
defined by 

S[yk{to)m:^yk{to)+ f F{yk{s))ds (3.55) 

Jta 

with yk(t) := T'^[xo]{t) and where T G C{C x M, C) is a succesive approximation 
operator analogous to the one defined in (|3.47l) . it can be seen that (|3.55p can be 
expressed in the form 

S[T''[xoKto)Kt) = T'^ixoKto) + f F{T^[xQ]{s))ds (3.56) 



from this equation it can be obtained that S[T'''[xo]{to)]{to) = T'°[a;o](to) wich 
implies that <5'[(-)(^o)](io) — it can be seen also that 

SiT^ixoWm') = SlSiT^lxoKtoWm') 

= T''[xo]{to)+ [ F{T''[xo]{s))ds + [ F{T''[xo]{s))ds 

Jto -It 

ft' 



= r'=[a;o](io)+ / F{T''[xo]{s))ds 

Jto 

= S[T'^[xo]itom') 
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Now if we take a time interval [to, r] to be partitioned in the form i„ to + nh, n = 
0, • • • ,N and where h = {t — to)/N the last relation presented above permits us 
to obtain a discrete expression for the time integration operator Sh[xn\(Tn) in the 
following way 

SMO)] S[T''[xo]{to)]in) - T'^IxoWo) + / i^(r'=[xo](s))ds (3.57) 



using this expression and without any lack of generality it can be observed that for 
any two m,n £ "Zq such that m > n we will have 

S„,oSn[xkm = Sn[xkm+ F{T''[xo]{s))ds (3.58) 

pto-\-7ih 

= T''[xo]{to)+ / F{T''[xo]{s))ds 

Ho+nh 



ta 



+ / F{T''[xo]{s))ds (3.59) 

•/ tQ-^-nh 

r''^N](to)+ / F{T^[xa]{s))ds (3.60) 



= S^+MQ)] (3.61) 

From (|3.57p and from (|3.6ip it can be seen that Sa[x\ ~ x and also that the set 
{S'„[T'^[-]] : n e Zg } defines a discrete semigroup with respect to composition. 
The action of discrete time integration operators over fixed elements of the approx- 
mate solution to the evolution equations can be computed using quadrature rules 
with accuracy order higher or equal to the one induced by succesive approximation 
operators. 

4. Approximation of Nonlinear Evolution Equations 

The nonlinear evolution equations studied in this work will be related in some 
suitable sense to elliptic operators not necesarily linear. 

4.1. Approximation of Semilinear Elliptic Problems. For a given semilinear 
problem of the form: 

A{u) = f{u),x£ G 



Bu = Ub,xedG ^"^-^^ 

where A{-) e jC{H"-{G)) is _ff " (G')-coercive exactly factorizable on £>/ {H^^ (G)) , one 
can obtain a particular representation of A in £/n^h{H"^{G)) using a particular 
projector compatible in some sense with the boundary conditions of m G H"^{G) in 
(|4.ip . Now, if we have that / is locally Holder continuous with exponent one ,i.e., 
/ G C"(i3, (0)), with Br{Q) a closed ball of radius r centered in 0, from here on we 
wiU define the set ^"(6; r) C H'^{G) by 

n"'{G;r) := {u e H"'{G) : u e B,.(0),Vx G G} (4.2) 

also if there is an < G M that depends on r, small enough, such that there 
exists c/ < oo with 

ll/(^) - fi^)\\H'^iG) < c/ \W - vWh^^^g) (4-3) 
when ||m — ^\\h"^(g) — then we can obtain the following. 
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Theorem 4.1. If for the semilinear boundary value problem 

A{u) = f{u),x e G 



Bu = Ub,xEdG 

we have that f G C"(-Br(0)) with a = 1, then we will have that the sequence defined 
by iiQ :— pj^^ ^uo,ma;+i := Gn.hf{uk) for a given uq G X{G) converges to a unique 
approximation of the solution to (j4.4l) determined by Uq . 

Proof. It suffices to show that under some particular conditions Qm.hf{') is a con- 
traction, it can be seen that for any u,v € X(G) such that Rg{u), Rg{v) C Br{uQ) 
for any x € G we have 

\\G.n,hf{u)-grn,hf{v)\\ < \\g.n,h\\ W f i^) ' f {y)\\ (4-5) 

< \\gm,h\\cf{r)\\u-v\\ (4.6) 

< 2M\\u^v\\ (4.7) 

m 

the last expression implies that r must be chosen such that c/(r) < m wich implies 
that f/„i, /!./(•) becomes a contraction, and by T I2. li the result follows. □ 

Theorem 4.2. // € X^^hiG) denotes the numerical solution to the particular 
representation of (|4.4p given by 

AmM'^h) ^ f{uh),x e G 
Bm,hUh = Pm^h'^b,x edG 

and if u := lim/j_j.o+ ^h, then we will have that 

\\U -Uk\\< Cuh-- + ^-2MMl{r \\AraM\\ + Mf)^Ji{G) (4.9) 

m — Cf 

where Mf := sup„g.^,i(Q.^^ |/(u)|, /i(G) ~ ||1|| and Vm is the approximation order 
of the particular projection. 

Proof. From R I2.il we have 

\\uH~Uk\\ < /%^'VT/\ H^i-^oll (4.10) 

^ /%^'"VT/\ l|grn../(^o)-ao|| (4.11) 

< /%^'VT/\ l|grn,^llll^»./»(ao)-/(ao)|| (4.12) 

< ("^^^)/^)' (||A„,,||||^o|| + ||/(uo)||) (4.13) 
(m - c/(r)) 

^ ^^^^")/7)' (P„,,||r||l||+A//l|l||) (4.14) 
(m - c/(r)) 

< [|+M;)MG) (4.15) 
(m - c/(r)) 
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also it can be seen that 

\\u~Uk\\ < \\u - Uh\\ + \\uh ~ Uk\\ (4.16) 

< \\u - Pm.huW + fi^'^^/"^)' (r WA^^^W + Mf)^i{G) (4.17) 

(to — Cf(r)) 

< cuh''- +i^^^^^^^^{r\\A^,h\\+Mf)fi{G) (4.18) 

(TO-c/(r)) 

from the last expression we get the desired result. □ 



Example 4.1. Given the semilinear boundary value problem: 



u 



= o,xea[o, 1] 



2 



(4.19) 



with A := d^+d'^, one can rewrite it in a particular representation of the form 

(02,1/32-^2, 1/3202, l/32)'«/i = X2, 1/32(1 + ufi) (4.20) 

where 02,1/32 S •!2^2,i/32(-f^o (G*)) is the particular representation of the exact factor 
0/^2,1/32 •= P2 1/32^^*2,1/32 o'T'f^ .^^2,1/32 the inner product matrix form relative 
to Hq(G). Taking Q2.1/Z2 to be defined by 

02,1/32 (02,1/32-^2, 1/32^2, l/32)-^2, 1/32 (4-21) 

and since a simple computation permits us to see that to — 2n'^ , we have 1 1^2, 1/32 1| 
{2Tr'^)~^ , on the other hand it can be seen that 

II/N-/(«)IIl2([o,i]2) = ||«'-^'|L2([o,i]2) (4.22) 

< (II«IIl2([o,i]=) + II"IIl^([o,i]^)) (4-23) 

< 2r||?/-z;||^.([o^j.) (4.24) 
now combining the results presented above we obtain 

||02,l/32/("l/32) ^ 02,l/32/(''^l/32)||^2([o.l]2) — l|02,l/32|| 

X ||/(wi/32) - /(Wl/32)||i2([0,i]2) 
< ^ ||ui/32 -l'l/32|L2([o,i]2)(4.25) 

SO we need to choose r < to ensure that Gm,hfi:) becomes a strict contraction, 
then by theorems T \4.1\ and T \4.2\ we will have that (|4.19p has a unique approximate 
solution ilk that satisfies the following estimate 

\\u -uk\\< c„(l/32)4 + i^Z^(ri^^(32)2 + 1 + r^). (4.26) 

4.2. Approximation of Dissipative Nonlinear Evolution Equations. In this 
section we will focus our attention in evolution equations of the form 

A{uit)) + f{u{t)) = u'{t) 

u{0)^uo,uoeH"'{G) ^ ' 
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where such that A{-) G j2/(iJ™(G)) will be an exactly factorizable operator on the 
particular C*-algebra £/{H"{G)) of the form 

[a'fV{v)a]v),x £ G 
Bv = vb,x e do 



A{v) 



with D{v) e a symmetric Tensor function, i.e. D{v)* — 

f G C°'=\n''{G; r)). 

In the following sections we will derive several useful estimates 



(4.28) 
D{v) and where 



4.2.1. Diffusion Tensor T>{u) IC. Using the above results, the first case that 
we will study is the one consisting of D{v) := /C, with K. constant in time and 
(SPD), in this cases and all the other we will study in this work we will consider 
that the boundary operator B G C{H'^{G)) allows A{-) G £/{H"{G)) to be exactly 
factorizable in .c/(iJ"(G)) and that / G C°'=^{n"{G;r)). 

If we take g G F{H^{Gr)) to be the right hand side of (|4?27| . i.e. g{v) := A{v) + 
f{v) and if we denote by A\I\ the form corresponding to A{-) when D{u) := I = 1 
then for any two u,v e dom{A{-)) n 'Hl{G;r) and since / G C='^{n"{G;r)) we 
will have 



\\giu) ~ g{v)\\ = \\A{u-v) + f{u)-f{v)\\ (4.29) 

\\Aiu-v)\\ + \\f{u)-fiv)\\ (4.30) 

\\a^lCa\\\\u-v\\ + \\f{u)~f{v)\\ (4.31) 

\\af\\lC\\\\u^v\\\\f{u)-f{v)\\ (4.32) 

||/C|| \\A[I]\\ w|| +c/ \\u-v\\ (4.33) 

^ \\IC\\\\A[I]\\+Cf)\\u-v\\ (4.34) 

hence g G C"^^(?^"(G; r)). In the following section we will generalize these results 
and present a useful technical theorem and its corresponding corollary. 



< 
< 
< 
< 



4.2.2. Diffusion Tensor V{-) G C%='^ 
the case when V{-) e G^=\nl{G;r), 



),C"^"). In this section we discuss 
in this particular cases we will also 



consider 



l£°°(G;r) 



also we will have that in these cases C*-identity is not 



satisfied in general, nevertheless these cases are also interesting in the analysis of 
nonlinear evolution equations, see example E 14.3[ 2?(m) is still (SPD) and is also 
bounded above by Mb := sup^g^fc^g.^) [[©(u)!! < oo. The conditions presented 
here imply that if w, w G dom{A{-)) n?^|(G;r) and if we represent A{u) by 
with A[u] := a^X'(u)a we will have 



\\A[u]\\ 
where fj.{A[I]) := \\a\\ 
\\A{u)~A{v) 



and also 
\\A[u 
\\A[u 
\\A[u 
\\A[u 
\\A[u 



< 



\a\\\\an\mu)\\<MMA[I]) 



A[v]v\\ 
A[u\v\\ 
i~v\\ + 



\\A[u]v ~ A[v]v\\ 
v\\ \\A[u]-A[v]\\ 
I \\u-v\\ + \\v\\ \\a'' {V{u) - V{v))a\\ 
I \\u-v\\ + \\v\\ \\a]\ \\a^\ \\V{u)-V{v)\ 
MDf^{A[I]) \\u ~ v\\ + CDrpi{A{T\) \\u - v\\ 
{MD^CDr)yi{A\X\)\\u-v\\ 



(4.35) 

(4.36) 
(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 
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from these we get the following estimate 

\\g{u)-g{v)\\ < \\A{u)-A{v)\\ + \\f{u)~f{v)\\ (4.43) 
< {MD+CDr)fi{A[I])\\u-v\\+Cf\\u~v\\ (4.44) 
= {{MD+CDr){fi{A[I])+Cf)\\u-v\\. (4.45) 

Hence g € C"^^ {C°° {G; r)). The last results can be sumarized in the following. 

Theorem 4.3. Nonlinear evolution equations of the form (|4.27p with diffusion 
tensor V{-) e (7^=1 (H^' (G; r), C"^") and with f e C"=^{H''{G;r)) have a unique 
local solution u{t) € Hg{G;r) in the time interval t G (—(5,(5) with 

" { rMz3M(^[I]) + Mf ' {Md + CDr)^iiA[I]) + c/ } " ^^'^^^ 
Proof From (|05)) it is clear that g e C'^=^{jC°°{G;r)) and since 

\\giv)\\ = \\Aiv)+fiv)\\ (4.47) 

< \\A{u)\\ + \\f{u)\\ (4.48) 

< rAf^M(^H) + 11/(^)11 (4.49) 

wich implies that 

Mg= sup \\g{v)\\<rMDKA[T])+Mf (4.50) 

tie-H^(G;r) 

now taking Cg :— {Md + cor) fJ,{A[X]) + Cf the result follows. □ 

Corollary 4.1. // Uh^n denotes the numerical solution of the particular repre- 
sentation of (j4.27l) after n iterations of the Picard method and if we take us '■= 
lim„_j.oo and if u is the solution to (|4.27p for t £ (—(5,(5) with d defined in the 
theorem T \4.3\ then we will have 

||„ - „,„.|| < + (((Mo + ...)MA[I|)+.,)fl" ^ 

1 - {{Md + CDr)ii{A[I\) + Cf)d 

Proof. If we take us P^njiU for some particular projector Pm,h G T^{T~l'\{G',r)) 
and from from corollary C I3.3I we will have 

\\U-Uh^n\\ < \\u - Us\\ + \\US - Uh^nW (4-52) 

< Cuh"'" 

, {{{MD+CDr)f,{A[I])+Cf)dr Mfrn^AfWA^^^ 

+ -j TTTT — \ X / .r-^n , V75{rii{A[I\) + Mf) 4.53) 

1 - {{Md + CDr)^i{A\I\) + c/)(5 

where Vm is the projection approximation order. □ 

4.2.3. Examples. In this section we will present some examples of implementation 
of the techniques presented before in the estimation of numerical solutions to the 
particular forms of some equations of mathematical physics. 

Example 4.2. Given a nonlinear Schrodinger Equation of the form 

—Au + jitpM = idtu, X £ G 

u = 0, a; € 9G (4.54) 
u{x, 0) = uq, uq G 'H{G; r) 
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where H(G; r) is a discretizable Hilbert space, if we take the particular representa- 
tion of (|4.54|) with respect to some Pm,h G 'P(H(G;r)) then we will obtain 



j u'^{t) = i[M^f^a*-^j^Mm,ham,h]'^h + iWh^^Uh 
it can be seen that 

and 

wich implies that 

m,h\\ \\Uh ^'i [^(Grr) 



fn.heinijh 



1 K 

< niMmM) \\am,h\\ < -J^ 



< 



K„ 



(4.55) 

(4.56) 
(4.57) 

(4.58) 
(4.59) 



^ \\uh ~ Vh\\^f^(j.^^^ 

then by T \3.& we will have that (j4.55p admits a unique local solution Um,h(t) for 
t e (-(5, 5) with 

1 1 



also by C\ 



t = min 
we will have 



Kjh^ + 3r^fi{G) ' Ka/h^ + r^ 



(4.60) 



(4.61) 



(4.62) 



\\u u^^f^wscun +^_^Ka/h^ + ir^tJi{G))5^ '''^ ' 
Example 4.3. Given the nonlinear diffusion equation 

dtu = dx{v?dxu) 
M(a:,0) = uo e n(G]r) 
uiO,t) = u{l,t) = 

where TilG^r) is a discretizable Hilbert space, one can obtain its particular repre- 
sentation with respect to some Pm,h £ T'{'H{G; r)) in the form 

i <it) = [M~,\dlj^Mrn.h{ul)d„,,h]uh 
1 ^pL.h^o 



(4.63) 



we have seen before that f{v) — G C" ^{T-L{G;r)) wich leads to the estimates 



\\AmJi[Uh_ 



< K{Mrn,h)\\dni,h\\ 

< r 



m.h 



r^llll 



/l2 

on the other hand we will have by ()4.45p 



K 

\\A,nAuk]uh - Aji./iK]w/i|| < 3''^-^ W^h - Vh\\ 



(4.64) 
(4.65) 
(4.66) 

(4.67) 



hence by T\4-.3\ exists a unique local solution Uh{t) of (|4.63p for t G (—(5,(5) with 



(5 = 



(4.68) 
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moreover by C \3.S\ we will also have that a solution to (I4.62|) will satisfy the estimate 

[ir'^KaS/h'^f r^KaS 



(4.69) 



1 - {3r^KaS/h^) h^ 

Example 4.4. For a vorticity transport problem in the vorticity- stream Junction 
approach that has the form 

dtu = JVv ■ Vu + ^Au 

" (4 70) 

Biu = Ub, X e dGB2v — vi,,x (1 dG 

where Re is a fixed Reynolds number and where JT^ G M^^^ is given by 

-1 



J = 



1 



(4.71) 



since we assume ^{Gir) is factorizable we can obtain a particular representation 
of (I4.70p with respect to some Pm.h S 'P{T-i{G]r)) given by 

taking -Bm,/ib], ^m,h G i^{G;r) to be defined in the form 

Bm,hM '■= Ja,n,hGmJiV ■ am,h (4.73) 



A. 



Re 



ft](am,/i) 



(4.74) 



we can rewrite (|4.72p in the form 

U'i^{t) = Bm,h[uh]uh + A„i,hUh 

and taking g{v) := Bm^h[v]v + Am.hV we can estimate Mg, first we will estimate 

\\BrnAv]\\ 



(4.75) 



\Bm.h\v] 



\J am.hQm.hV ■ am,/i|| 



< ||0m,/i|| \a.mJi\ 
KaT 



\v\\ 



< 



m 



so M„ will satisfy 



^ ' Re 

on the other hand we will have that 



\\g{u) - .g(t;)|| < \\B^Au]u - B„^Au]v\\ + \\B,nAu]v - Bmji[v]v\\ 
+ \\A„,.h{u-v)\\ 

= lj(^ + ^<Mrn,H))\\u-v\\ 

h'^ m Re 

so it can be seen that g G C"^^('HJ^ i^{G;r)) and we can use T \3.6\ and C\ 
perform similar estimates to the performed in examples E\^.S\ and E\4-.3\ 



(4.76) 
(4.77) 

(4.78) 
(4.79) 

(4.80) 
(4.81) 

(4.82) 
to 
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